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I. INTRODUCTION 

The lack of analytic solutions for the most centrally 
symmetric potentials calls for developing method approx- 
imating the spectra. It becomes extremely actually at 
present as besides known problems in atomic physics, 
theoretical models of nuclei [l[ or quarkonium 0, H, 0] 
there arises a number of new objects (metallic clusters 

[f| etc) for which some analogies of the periodic system 
of shells may be constructed @, d, 0] . They may differ 
not only by the nature of selfconsistent field but even by 
their dimensionality d. A lot of endeavours have been put 
to obtain rigorous results, see e.g. 0, H, H [2(| HH, [22j], 
but these theorems only give some inequalities for special 
forms of potentials and only for d = 3. 

A frequent supposition is that the energy values de- 
pend on some linear combinations of the radial n r and 
orbital I quantum numbers, i.e. E(n r , I) — E(an r + (31). 
The Madelung-Kletchkovsky rule predicts appearing of 
new shells (n r ,l) in the Periodic system of the elements 
with increasing n r + 2l 0,[lJlKfor the metallic clusters it 
is expected E = E(3n r +l) [5|,l3K a similar quantum num- 
ber for nuclei was proposed in |l[ by using some classical 
analogy. 

However, such dependence of the exact spectra on 
some linear combination is known only for the Coulomb 
and oscillator potentials: 



V c (r) 



Z 



E c (n r ,l) = -C c {v + \y 



(1) 



V osc (r)=br 2 , E osc (n r , I) = C osc [y + A/2) , (2) 



i/ = n r + i, \(d) = l + ^^. (3) 
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Two remarkable facts are known (at least for d = 3) for 
our reference potentials ([I]) and ([2]) and only for them. 
First of all, the usual WKB condition 



~ / \j2(E-V{r))-^dr = V 



(4) 



leads to the exact spectra (dJ, @. For d = 3 the 
term A 2 /r 2 is known as the centrifugal potential with 
the Langer correction 1(1 + 1) — » (I + 1/2) 2 [17]. Besides, 
for the oscillator 



2(E-V(r)) 



A 2 



-dr = 



- f ^2(E-V(r))dr 



(5) 



In all cases the integration goes between corresponding 
turning points and we put m = h = 1. Combining ^ 
and (0, we obtain 



i J ^2(E-V(r))dr = v + ±, 



(6) 



so that in the right hand of ([6]) we see the same linear 
combination of v and A as in @ , a similar situation with 
v + A we obtain for |1| . 

In the present paper we get an effective quantum num- 
ber 



T (n r , I) = v + 4>\ 



d-2 



(7) 



for any centrally symmetric potentials and any dimen- 
sionality of the problem by means of some generalization 
and new universal variations of the methods used ear- 
lier El Git in order t0 obtain T for d = 3. The 
coefficient <f> is determined as some definite combination 
of the functionals Nd[E;V], which represent asymptotic 
estimates for the total number of the bound states in a 
given potential V(r) with the energies not exceeding E. 
These estimates are obtained without using any WKB 
method. 
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This quantum number T ||7J determines the order of 
the bound states: E(T) > E(T) if T > V . For our 
reference cases (JTJ and ((2]) T coincides with v + A and 1/ + 
A/2 correspondingly. In a general case T approximates 
very good the real situation, but is not formally the exact 
quantum number. 

Note that the principal quantum number n = n r + 1 + 1 
determines the spectrum only for the Coulomb potential 
if d = 3. 

II. THE SCHRODINGER EQUATION AND ITS 
CONFORMAL TRANSFORMATION 



We use the Schrodinger equation in the form (to = h = 

1) 

A d <S> + P 2 $ = 

/ d 2 d-ld\ x A d _i(fi)$ , n 

= I AT + IT $ + 2 + P $ = °. 

\ or r or J r z 

P 2 = 2(E- V(r)) , (8) 

where A £ j_i(J7) denotes the Laplace operator on the unit 
sphere S , P is the energy value and V(r) is the po- 
tential. The equation {SJ corresponds to the following 
metric: 

ds 2 = dr 2 + r 2 dn 2 d ^ (9) 

with < r < oo and f2<j_i being coordinates on S 

In order to eliminate the first derivative in ([5]) and 
the singularity at r = we represent the Schrodinger 
equation in a conformal metric ds 2 with a new variable 
p = In r so that 

ds 2 = r 2 d? 2 = e 2 " (dp 2 + dO^_i) . 

Using [9, §28] we get in the new metric instead of ([5]): 

( ?-^ + A d ^ 1 (n)y + K(d)y + e 2 PP 2 «f = 0, (10) 
dp z 

K(d) = -{d - 2) 2 /4; * = e^ p $, -oo < p < oo. 

(11) 

Certainly it is possible to prove dTUJ) by a simple substi- 
tution $ (jTTJ) into ||5J|; for d = 3 it is the Langer trans- 
formation introduced earlier ad hoc. 
Taking into account eigenvalues 

A d _!(f7)y = -L 2 Y; L 2 = l(l + d-2), 1 = 0,1,2,... 

(12) 

and the term K{d) in ([TT]), for * = ip(p)Y(Sl) we obtain 
ip" + Wtp - \ 2 tp = 0, A = / + ^-^, (13) 

W(E, p) = r 2 P 2 (r) = 2e 2 " (J5 - 1/(6")) . (14) 



Usual condition $(r = 0) < oo leads to i(>(p) — >• if 
p — > — oo. Exact spectra of © and PH|) . p3]) must be 
identical as $ and \& are only distinguished by a positive 
factor e 2p . 

We have seen that the usual "automatic" replacement 
1(1+1) —y A 2 actually means that we work in a new special 
conformal curved space; its curvature is Ke~ 2p with K 
from (flT]). Instead of the topology R d = S^ 1 x (0,oo) 
we get S ,d_1 x (—oo, oo) in the conformal space. 

Coordinates (p, f2) are similar to the Cartesian ones 
in the maximum possible measure: p is a harmonic co- 
ordinate, there is a field of the parallel vectors and all 
sections p = const are identical ones [l8|. That is why 
the leading WKB approximation in the conformal space 
gives the best possible result (while exact spectra must 
be identical in two metrics ds 2 and ds 2 ). 

It is known [ll| that the number of eigenstates with 
the same value of I is 

fl(0 = n ( 7_")M n, = (i + d-i)(*+d-2)...(i + i) 

(15) 

and 11; = if I < 0. Substituting in (JT5J) Aflwe get the 
leading term D: 

D(X) = D(X) + r(A, d); D = — — , d > 2, (16) 

where r(A, d) = if d < 4 and r is of order A d ~ 4 if d > 4. 

The leading term of I? as a function of I is simply D(l), 
but r(l, d) is of order l d ~ 1 . For example, if d — 4 we have 
D{\) = D(X) = A 2 = I 2 + 21 + 1 = I)(l) + 21 + 1. Thus 
in D(X) there is no term of order A d ~ 3 unlike D(l), so 
that namely A is the most suitable variable with smallest 
distinction between D and the leading term D. 

The WKB quantization condition for (fT3")) is 

I(E, A) = i J ^W(E,p)-\ 2 dp = n r + X - = v. (17) 

Obviously (fT7| is identical to (H]) when we return to pre- 
vious variables r, P 2 ©. 



III. LINEAR IN A APPROXIMATION OF THE 
WKB INTEGRAL 

As we have already seen, I(E, A) is linear in A for the 
oscillator and Coulomb potentials. In general case we 
can write 

N X (E) ee I(E, 0) = I(E, \) + (f>\ + q(E, A), (18) 

where q denotes all non- linear corrections. For above 
cases (P), (|2|) q ee and 0=1 and 0=1/2 respectively. 

For determining cj) we use hereafter the known esti- 
mates Nd for the total number of bound states in d di- 
mensional problem with the energies not exceeding E 
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[12j : being expressed in our variables W, p they look as 



N,, 



M d (E); 



n(d- 2)! 

M d {E) = [w d ^(E,p)dp, d>2, (19) 



where B(y, z) is the beta function. Let's show how N d 
(fT9)) may be obtained using (jTTJ) and (|16p. At a fixed A 
([5]) the equation (JTTJ) determines energies of the bound 
states. Evidently the maximum value of v corresponds 
to the maximum energy, so that I(E, A) is equal to the 
total number of bound states with a given A and values 
of energies not exceeding E (as usual in WKB we neglect 
the difference between n r and v in this case). Now we 
have to take into account the degeneration of states (fT6|) . 
Using the universal form of the first term D, we multiply 
I(E, A) (fT7|) by D and integrate with respect to A over 
all possible domain < A < A, 

A 2 (E) = maxW{E,p). 
p 

We suppose that W(E, p) has sole maximum at any E. 
In intermediate calculations we treat A as a continuous 
variable. Changing the order of integration we really 
obtain 



A{E) 



D(X)I(E,X)dX = N d (E) 



(20) 



with Nd from (| 1 Qjl . Thus namely the leading term D but 
not D must be used in WKB methods. 

Now we return to (|18[) and intend to choose value of <p 
so that q(E, X) averaged over all bound states becomes 
zero. Multiplying both sides of (JT5]) by D and integrating 
with respect to A we obtain 



Nd 
A d 



2N X 



20 



A(d-1)\ d{d-2)\ 



(21) 



Obviously N d is proportional to A d so that the value of 
4> does not depend on A. It is also invariant under the 
transformation r — > etr, i.e. p — > p + In a, a > 0. 

In order to simplify (|2ip we introduce relative param- 

(c) 

eters Xd comparing Nd for a given potential with N d 
for the above reference case: the Coulomb potential with 

(c) 

the same value of A as in Nd- Calculating N d we get 



N d 



Xd 



N 



d\N d 
2A d 



d>2: 



Xi 



* , 22 , 



We have used here some identities between the beta func- 
tions and their products. Substituting l|22[) in (f2Tjl and 
using (fTO)) we obtain finally 



= Xi + 



Xi ~ Xd 
d-1 ' 



,4 d m d ' 



m d 



5(d/2,l/2). 

(23) 



In general case we have = 4>(d,E). For our reference 
cases CD), © </> = Xd = const does not depend on d 
as well as on E and coincides with the discussed above 
= 1, 1/2. Thus we have received the desired effective 
quantum number T ([7]) with of a ([23]) . All the pa- 
rameters Xd, 4> are some functionals of a given potential: 
(j) = <j)[V] and so on as well as some functions of E and d. 
Hereafter we indicate only actual in each case arguments. 

The quantization condition (fT7f with respect to l|18[). 
(f2"Tj) takes the following form: 

1 



A(E)xi(E) = -j ^2{E-V{r))dr = T (n r , I; cf>(E)) . 

Functions Xd(E) may be treated as a special non-linear 
transform of a given potential V(r). This transformation 
is the most adequate one for our method. 

Note that 4> and Xd are invariant under the transfor- 
mation W -> BW with B > 0. 

In what follows we treated (|23|) as the basic form for 
cj). Nevertheless it is useful to find another approximation 
for <p and to compare them in order to evaluate their ac- 
curacy. Suppose we construct a multiplicative expression 
for (f> 

cjy m {E,d)=CN2N^ 

using only integral estimates N d and excluding the 
above parameter A. Since N d oc A d we must put 
ad + (3 = 0. Substituting Xd instead of N d we get 

(j) m = G (xiXd 1/C( ) • Demanding that <p m = Xd = 4> 
for our reference cases we have finally 



<>,„(£.</) = ( — 

\Xd 



(25) 



The same value of <p m we have received earlier using a 
duality between certain pairs of power-law potentials 14]. 
As we shall see, Xd is a smooth monotonic function of d 
for all interesting model potentials. It's conveniently to 
write 



Xd = Xoo exp/ d , 



(26) 



introducing an auxiliary function f d ; then the ratio of < 
and d> 



2(d-iy 



(fi-fd) +0(f). (27) 



Numerical calculations really give R > 1 but very close 
to 1, so that R — 1 < 0.01 and in most cases < 0.002 in 
a wide interval of E for a wide set of studied potentials 
and for all d, see Table |TJ 

This closeness of two values of cj) calculated by quite 
different methods confirms the objective character of <p 
and of the effective quantum number T itself. We shall 
return to this closeness in the section IVII1 

Note that we may treat d in M d ,Xd and so on as a 
continuous variable in intermediate calculations. 
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IV. NON-LINEAR TRANSFORM X d FOR THE 
POWER-LAW POTENTIALS 

In the present section we study power-law potentials 

V{r) = br^, bu>0 (28) 

with — 2 < \l < oo. For them functions Xd (22) and <j> 
(|23p are monotonic and not depending on E 

' 2 + ii \^i 2 d ' 2 Bjd/j^d/2 + T) 



f 2 + u \^^ 2 d ' 2 B{d/u,d/2 + l 
XdW \ 2 ) ^/ 2 +! B(d/2,l/2) 



(/i > 0) 



(29) 



Xd(^) = 



3 — It- 



V 2-| M |; | M |«»/a+i S(d/2,l/2) 

(-2 < /i < 0). 

In the limiting case d — > co we obtain from (|28| , ([29| for 
all /X > -2 

Xoo = -t=L=. (30) 

Naturally (|2"9"|) does not depend on d if u = 2, /i = — l. 
Values of Xd an d (J2S1) are shown in the Table l. 

For the power-law potentials the following expansion/^ 
([26]) with Xoo (TO is valid 



k=l 



(/j + 4) 2 
12(^ + 2 

1 
360 



3 

' 4' 
8/x 3 



(31) 



(M + 2) 



6 2 = 6 4 = 0. (32) 



The values of b\ and 63 are equal to zero for the oscillator 
and Coulomb potentials and are small in their vicinity 
with 1 63 1 <C 1 61 1, b\ 63 < 0. Correspondingly the ratio 
R (|27p is very close to one. At any fixed d the value of 
</>(d, n) (|23[) with Xd (|2T)|) decreases monotonically with 



increasing fi. 



V. THE LEVELS ORDERING IN 
ACCORDANCE WITH T AND EXACT RESULTS 

Since the levels ordering is a very important property 
of many-particle systems many authors have tried to ob- 
tain exact theorems for several special forms of poten- 
tials and only for d = 3. On the other hand, the effective 
quantum number T ([7]) together with the condition (I24p 
immediately leads to the universal for all potentials V(r) 
ordering 



sgn (E(T) - E(T')) = sgn (T - T') 



(33) 



since the left side of the condition (|24|) is a monotonous 
function of E. In the present section we show that our 
expression ((7J, (j2"3"|) for T is exact enough to reproduce 
many known or expected results. 

It is interesting to compare (I33[) with known results. 
The first example connects two differential operators 
(d = 3) 

= r-^(rV)=r^f(K(r) + l), 
dr* dr 



Y 2 [V] = 



dr 



n(r) = 1 + 



1 dV 

r dr 
(d 2 V/dr 



1 dV , , , 

(34) 



d^/dr 

with two energy differences correspondingly: 

Di = % + l,I)-%,! + l), 

L> 2 = E(n r ,l)- E(n r -1,1 + 2). (35) 

We have introduced an auxiliary function n(r); for power- 
law potentials fc(r) = /i. 

The theorem proved in [2(| [23| states in our notation 
that 



sgnYi — sgnDi, 
sgnY2 = sgnD2 



(36) 



if sgnYk = const for < r < co. Suppose our T- method 
is exact. Then, taking into account that for monotonic 
attractive potentials dV/dr > in (f3"4"| and using (|3"3"j) . 
we can write instead of (1361) 



sgnYi[V] = sgn(n(r) + l) 

= sgn(T(n r + l,l)—T(n r ,l + l)), 

sgnY 2 [V] = sgn(n(r)-2) (37) 

= sgn (T(n r , I) — T(n r — 1, 1 + 2)) . 

For power-law potentials <f> — c/)[u] does not depend on 
E, k — u and substituting T ([7]) into (|3"T|) we get two 
equivalent to (|36[) equations 

(38a) 
(38b) 



sgn(/i + 1) = - 4>), 

sgn (/! — 2) = sgri(l — 20). 



The equality (|38aj) really takes place for all potentials 
as <fi (1231) with x<i (12^1) is a monotonic function of 
fx, see examples in Table U ,0 = 1/2 for k = /1 = 2 and 
= 1 for k = fi = — 1. 

Another exact statement is the following one 0]: for 
power- law potentials ^ d 2 E(0,l)/dl 2 > if a > 2. 
From our condition (|24[) we get immediately 

E(n T ,l) = C v T^ =C ll (v + <f>\)&' , C^>0, u>0, 



sgn- 



d 2 E 



or 2 



sgn(u - 2) 
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so that the abovementioned inequality really fulfills with 
our T for all n r . 

As the third example we study a family of potentials 
for quarkonium systems Q 



V q (a,5,r) = £(-- + (l-a> 

< a < 1, S > 0, £ > 0. (39) 

The following ordering was discussed for V q 

£(0,1) < £(1,0) < £7(0,2). (40) 

For this family 9V/9r > 0, dfc/dr > 0, 



1 < «(r) 
Q(r) 



-i + Qm 2 

l + Q/x 

1 - 



< k(oo) = (5, 



(41) 



for any r > 0. Thus the above theorem (f3T>|) holds and 
predicts the ordering (j4"0|) , but only for r5 < 2 (if 5 > 
2, sgn(n — 2) changes its sign at large r). 

Now we study how our T-method works in this situa- 
tion. Since we act within the WKB frame it is natural a 
"weak" 

Conjecture 1. The levels ordering satisfies (|37p if sgn{n — 
2) = const only in the classically accessible domain: 



r < r t , 



V(r t ) = E. 



(42) 



(we suppose dV/dr > for our potentials). This con- 
jecture allows us to have Ac(rt) and correspondingly the 
previous levels ordering as for /i < 2 even 6 > 2. It is 
easy to calculate that for (f3"9"| dn(r t )/da < 0, so that 
we can have n(r t ) < 2 for small values of energy or for 
middle values of a but n(r t ) > 2 and the inverse ordering 
£(1,0) > £(0,2) if £(1 - a) is great enough. Rigorous 
results confirm it 0]. 

Besides, we have d(n m )/dE > 0. Note that for a <C 
1 the Coulomb term is negligible in (f39|) even for the 
deepest levels as compared with £(1 — a)r s . 

Moreover, it seems to be reasonable a "strong" 
Conjecture 2. For any potential with a smooth dn/dr the 
value of <p can be well approximated as cf>\p] for power-law 
potential (|28p with /i = n(r rn ). Here r m depends both 
on the value of the energy and the parameters of a given 
potential, and r m is defined as W(r m ) = max. 

Really, this conjecture is exact for Xoo, compare (|30|) 
and (|4"9")l . Each \d and thus </>, as it follows from (|3~lj) and 
(|47l) , includes functions depending on some derivatives 
lV^(r m ). We can express W^(r m ), a k (|4"T)) and corre- 
spondingly 6^ (f3Tj) by means of ft(r) (|34[) . For example, 
the main terms take the following form 



Xd = Xoo[^m] 1 + 



dil 



with fj, 
from 



K(r„ 



Xoo 

a new function 



accordingly to (O or (gH]), 

depending on n(r m ) in 



such a way that 6° = if «;(r) = const and 



fix 



dn{r m ) 
dr 



Using (|49l ) we can write for d = 3 with the same accuracy 

(f> = 4>{[J-m) + 4> add , 



idd 



46f 



0(1(1 



3V^m + 2 



with one of the previous expressions for 4>. In the most 

6. So linear in K r 



iadd 



or even all real cases we have 
and neglecting k^, n rr approximation is 



Using (j4Tj) we get 

dn 
dx 



16 + K(r m ) dn(r m ) 



24 (K(r m ) + 2) 



1 



rife 



fir Q(r m ) 



dr 



< 1 



at large Q. The condition Q ^> 1 fulfills at large £ 
and/or large B(l — a). Thus " non-adiabatic" correction 
< h, cj) add < 0. 
All the investigated potentials have \d, 4> as a smooth 
function both of £ and d as well as of other parameters. 
It seems to be correct for all physically well-founded po- 
tentials. 

Then the coincidence of the two levels 

T(0, 2) - T(l, 0) = 1-20 

is possible if 8 > 2 and £ or B(l — a) is large enough. 
So the line 1 — 2<fi = on the plane (6, B(l — a)) divides 
two domains with opposite levels ordering. Note that 
T m < ft] for the above case 



n < r t . 



2 + 6 



VI. SCREENED COULOMB POTENTIALS 

In the present section we study with the help of our 
new method another actual and interesting class of cen- 
trally symmetric potentials 

V (r)=-?M; ff (0) = 1, g>0, ^<0. (43) 
r dr 

The Thomas- Fermi potential VtfW) of the selfconsistent 
field in the many-electron atoms [lj| also belongs to the 
type (|43[) . All such model atomic potentials must obey 
the inequality g" > 0. It follows from the Poisson equa- 
tion 



AU 



Yi[U] 



-47rp 
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TABLE I: Values of \d and <j> J23]) for several potentials: (|4"5kb.c). (|45H) Vtf(t-) and (J28J) , as well as m (3) ([25]) for these 
potentials. 





V(r) 




X3 


X2 


Xi = 0(oo) 


0(3) 


m 


0m(3) 


E = 




1.414 


1.376 


1.359 


1.316 


1.286 


1.273 


1.286 






2 


2 


2 


2 


2 


2 


2 




(HSb) 


1.826 


1.803 


1.793 


1.769 


1.752 


1.745 


1.752 






1.89 


1.87 


1.84 


1.78 


1.74 


1.72 


1.75 


£ -» -oo 


(|45b.b.c,d) 


1 


1 


1 


1 


1 


1 


1 




ii = -\ 


1 


1 


1 


1 


1 


1 


1 




V -> 


0.707 


0.688 


0.680 


0.658 


0.644 


0.636 


0.643 


V = &r" 


/x = l 


0.577 


0.568 


0.563 


0.551 


0.543 


0.539 


0.544 




^ = 2 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 


any E 


/i = 3 


0.447 


0.457 


0.461 


0.469 


0.475 


0.477 


0.476 




/l — > oo 





0.212 


0.250 


0.318 


0.371 


0.386 


0.390 



for the electrostatic potential U, so that V = — \e\U and 
the electron charge density p < 0. Obviously sgnYi[V] = 
—sgnYi[U] = — 1, correspondingly we get n < — 1 from 
(f3~4"| . On the other hand, for (|4"3")) we obtain 



k= -1 



g"r 2 
g'r- g 



= -1- 



<A 2 



1.91 



(44) 



with respect to the inequalities (|4"3")) , so that actually <?" > 
0. 

There are also important potentials with 



9{r) 
9( r ) 



(1 + r) 2 ' 
1 

(1 + r) 2 - 5 



(45a) 
(45b) 

(45c) 



besides (gSbl) Vrf (r). 

For all these potentials the values Xd and 0(d) depend 
on the value of the energy E. For the deepest levels where 
only a small domain r < rt — > is classically accessible in 
(jMj) (V(r t ) = £J), we have k -» -1, and Xd -» 0(d) -> 1. 
Of course in our quantum problem the deepest level has 
some small rt ^ as well as E ^ — oo, so that formally 
we have a very small distinction from this limiting values 
of Xd, 4>{d). 

In the opposite extreme case E — ► — correspond- 
ing values Xrf and 0(d) for (j4"5a|) , (|45b|) . (|43c|) and ijlSH) 
Vtf{t) are shown in the Table |TJ These values asymptot- 
ically coincide in the depth of the potential well where is 
no screening as the turning point rt —* for these V(r). 
In the same Table we have placed Xd and 4>(d) for some 
power-law potentials (|28p , including /i — > oo (the rectan- 
gular potential well) as well as multiplicative <j) m (3) from 

(ESI. 

Let's demonstrate two ways to use T ([7]) taking as an 
example the model potential ([43b ) which is a very good 
approximation of the real selfconsistent atomic potential 
[15| with Z being the nuclear charge. If we fix E = 
then (|24p indicates the order in which new bound states 
appears with increasing Z as well as corresponding values 



of Z . On the contrary at Z = const we get values and 
thus the level succession of all the bound states in a given 
atom with a fixed Z . It can be easily seen from our ta- 
ble and (|2"4"1) that shallow levels (E w 0) are governed by 
n r + 1.751 (of course d — 3) but deepest levels by n r +l i.e. 
are Coulomb-like, with an intermediate behavior for mid- 
dle levels. The introduced T-ordering formalizes the well 
known quantitative picture and in particular explains the 
Periodic system of the elements, see the next section. 

Our method is valid and simple not only for potentials 
with non-trivial analytic form as (|45c|) but also for poten- 
tials given numerically. So for the Thomas- Fermi atomic 
potential we obtain 0(3) very close to the corresponding 
value for (145cD. 



VII. UNIVERSAL DIAGRAM 

The regular filling of shells in a centrally symmetric 
system with an arbitrary dimensionality and nature of 
the selfconsistent field is clearly described by the follow- 
ing diagram Fig. [T] Here each line represents T(n r , I, cf>) 
as the linear function of at fixed (n r ,l). Crossing of 
two lines marks the values of at which the order of the 
level succession changes. Two different types of problems 
may be treated with the help of this diagram. 

If E = the value of is invariant under the transfor- 
mations V —>■ cV, r — * cir, this value does not depend on 
Z (even if g = g (rf(Z)); that's not the case if E ^ 0). 
With increasing strength of the potential (i.e. Z in ([4^]) 
) new shells (n r , I) with E = appear in the order of in- 
creasing T(n r , I) at a given value of <j>(E = 0), i.e. along 
the vertical line = const. Each shell can contain D(d) 
([Tn]) states; only this number D depends on d. When the 
shell (n r ,l) is filled in, i.e. all states are occupied with 
particles, begins filling of the next shell (n' r ,l') with the 
value T' >T nearest to T. 

It is known that the structure of the Periodic 
system of the elements corresponds to the def- 
inite order of the atomic shells filling with in- 
creasing Z 0, [HI]. The actual shells order is: 
Is, 2s, 2p, 3s, 3p, 4s, 3d, 4p, 5s, id, 5p, 6s, 4/ i. e. 
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FIG. 1: The effective quantum number T © versus 4> p3jl 
for the bound states (n r ,l) with n r and I being the radial and 
orbital quantum number respectively. Lines with I > 3 are 
omitted for the sake of simplicity. 



(0, 0), (1, 0), (0, 1) etc. It takes place if and only if 

5/3 < 0(3) < 2 (46) 



value 0(3) « 0.37 for all levels we obtain for the rectan- 
gular potential well (|28p with ji — ► oo. 

We can also add that the T-ordering or immediately 
Fig. [1] exactly reproduces the levels ordering calculated 
for (|28|) with d — 3 and several values of e. g. for 
H = 0.1, « 0.63 [111 and /x oo, » 0.37. 

Another situation arises if we treat 0(2?) and the left 
side of the condition (|2"1)) .A(2?)xi(2£) at a fixed £ as a 
point (0(2?), A(E)xi{Ej) on the plane of the diagram. 
Changing E we receive a curve, the form of which de- 
pends on a given potential V(r). 

At an arbitrary point we can't find integer (n r , I) satis- 
fying ([24)1 . Such integer (n r , I) only exist for these points 
where one of the lines T(n r , I, 0) crosses the curve. Thus 
these distinguished points indicate the actual bound 
states and indirectly their energies. The curve as the 
whole shows clearly how the levels ordering changes with 
changing E or T (remind that dT/dE > 0). 

As an example we indicate on the diagram such 
curves for the Yukawa potential V(r) = — 50e~ r /r 
(the right curve) and the quarkonium potential V(r) = 
3 (— 1/r + r) (the left curve) . 



VIII. ASYMPTOTICAL BEHAVIOR OF Xp AND 
AND A NON-LINEAR QUANTIZATION 
CONDITION 



As we have demonstrated in previous sections, our 
method with the linear in A approximation for T is suf- 
ficiently exact. Meanwhile, it is possible to take into 
account non-linear corrections. In the present section 
we introduce a non-linear form for our effective quantum 
number T and demonstrate its new possibilities. Let's 
use the fact that the integral (fll?)) has the original form 
for the asymptotic Laplace expansion: 



Md 



dxW d/2 (x) 



A d+i 



4tt 



\WM\d 



V d d? 



(47) 



where a& are known functions of the derivatives 
taken at the maximum point r = r m , A 2 — W(r m ). 
Expanding in the denominator of (|23p similarly to 
(l47|) we obtain an asymptotic expansion of \d (|26p . (I3ip . 
where 



(48) 



as it follows from Fig. [TJ for the first time this range 
was found out in another equivalent form of the atomic 
potential asymptotics [lj, [l5| • The most probable value 
(j> w 1.75 corresponds to the Thomas-Fermi potential as 
well as to (|45cp and satisfies pp) . 

The usually supposed shells ordering in metallic clus- 
ters @ means in our notation that <fi = 1/3. The close 



A special case of such expansion is (|32|) . Of course we can 
also calculate Xoo as lim Xd for d — > oo (even for nonana- 
lytic potentials, when l|48p may be incorrect). It is easy 
to prove numerically that all interesting model potentials 
including those represented in the Table U have b\ <C 1, 
I&2I I fox I and so on. E. g. for power-law potentials |28p 
we have accordingly to §2$ h(0) = -0.08, h (4) = 0.139 
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and only in an irreal case n ss 15 we get b\ (//) « 1. Thus 
the asymptotic regime is already reached for d > 1 (actu- 
ally for d > 0.5). That's why Xd are smooth monotonic 
functions for d > 1. 

Using At(r) ([34]) we can get an interesting expression 
equivalent to ([48]) : 



Xoc 



«(r m ) + 2' 



(49) 



For power-law potentials k = u and we return to (|30[) . 

Neglecting in ([2"6"]l . ([3T]) all terms with d > 2 we 
obtain a simple approximation 



~ ( as ) 
Xd ~ X d =Xo 



Xi ~ Xo 



Substituting ([501) 

in (|23[) . we also have 



</»(d)«0( as )(d)=Xi + 



Xi - Xc 



(50) 



(51) 



At last comparing ([50]) and ([5T]) we obtain another 
approximate expression 

d 



d+l 



(52) 



We have calculated ratios s = 4>^ as ^ (d) / <f>(d) and 

Xd S V < / ) ( c where 0(d) is the basic form ([23]) for vari- 
ous potentials and for a wide interval of E. It turns out 
that even for d = 3 both s and w are close to unity: 
\s — 1|, | io — 1| < 0.02 and in most cases < 0.01. 

As we have already said, the closest approximation to 
the basic additive form cj> is the multiplicative form </> m 
([25]) : for their ratio R = <fi/<f> m we have R — 1 < 0.01 and 
in most cases < 0.002, see Table Q] 

Each of these approximations may be preferable in 
some suitable situation. 

A simple universal form xi| ([511)1 allows us to get an 
universal non-linear approximation 



I(X) = 1(0) - F(X) 



(53) 



for I(E, A) (jXTJ) . Acting as in Sec. IIIII we obtain instead 
of 



A d 



F(X)X dX 



d-2j, _ dxi-Xd 



d-1 



(54) 



If we assume F = (j>X in (J5J|, we return to <f> ([2"3"]l . But 
this time we use Xd (|50p and thus come to a simple 
Mcllinc transform for F(A)A _1 : 



A d 



F(X)X dX = 2± + 



Xi , Xi - Xo 



d 2 ' 



F(X) = X iA+(xi-Xco)Aln(A/A). (55) 

Correspondingly we obtain a non-linear in A quantization 
condition 

A(E) X i(E)=T non (n r ,l), l) = v + F{X). 

(56) 



Both ([55]) and (| 56 j) become incorrect if A/A <C 1, but 
actually if d > 3 we have a finite A/A even for / = 0, see 
(| 1 3[> : in real systems usually A A maa; < 5. 

Let's reproduce a delicate distinction obtained for the 
power-law potentials (1281) b y means of the ^-expansion 
of the Rcgge trajectories 23]. Using our notation we can 
rewrite the final result 23] in a simple form for [i > — 1 



sgn (E a (0, 1 + 1) - 2E a (0, 1) + £^(0, 1 - 1)) 
2^i 



sgn(2-fj,), a 



H+l' 



(57) 



Since E cx for ||28j) the left side of ([57]) is linear in T; 
by substituting T non ([55]) in ([57[) . we obtain 

S5« (Xoo - Xi) ' sgn A = sgn (xoo - Xi) 

= 8 gn(2-ii), (58) 
A = (A + 1) ln(A + 1) + (A - 1) ln(A - 1) - 2A In A. 

We have used here the fact that d 2 (x In x)/dx 2 > 0, so 
that A > 0, sgn A = 1. But the last equality ([58]) really 
fulfills for our Xd, see e. g. Table [H Should be stressed 
that in any linear approximation the left side of Eq. (1571) 
is equal to zero, and not to one as for A, so that onl 
non-linear approximation confirms the strong result 



IX. CONCLUSION 



only 
23]. 



Thus we have constructed and calculated the effective 
quantum number T ([Jj which determines with very high 
accuracy appearing and ordering of the bound states in 
any centrally symmetric potential. It should be stressed 
that different potentials may have very close or coinciding 
values of <j> and T and hence identical levels ordering (see 
e. g. cases b) and c) in the Table [J). 

Some partial success of using the Thomas- Fermi poten- 
tial to explain the Periodic system does not mean that 
this potential is the genuine or the best one. The point is 
that its 4> is situated not too close to the limiting points of 
the interval (|4"6")l so that various corrections not taken into 
account in the Thomas-Fermi approach though change <fi 
but remain it inside this interval. Hence the levels order- 
ing is really the same as for the Thomas-Fermi potential. 

The effective number T actually replaces the principal 
quantum number n = n r + 1 + 1 for all potentials besides 
the Coulomb one with d = 3. 

Using T we immediately reproduce many results re- 
ceived for the levels ordering by means of special theo- 
rems and numerical calculations. Moreover, the quan- 
tization condition ([24]) with T can also be used for de- 
termining spectra. The accuracy of the energy values 
calculated for V(r) — r is 0.3 -j- 0.5% and even in the 
worst case of the non-analytic potential well (/x — ► oo) 
errors do not exceed 3-^5% 16]. For two potentials with 
the same value of <p we usually obtain different A and xi , 
so that their energies of the bound states do not coincide 
unlike the levels ordering. 
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A wide variety of the physically interesting poten- sal asymptotic behavior if we use special adequate trans- 
tials, quite different as transforms of r, have an univcr- forms \d and 4>. 
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